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ABSTRACT 


The Fram Strait is a region of complex circulation and intense dynamical activity. 
Its general circulation is largely influenced by the East Greenland Current (EGC), the 
associated Marginal Ice Zone jet, the West Spitzbergen Current (WSC) and the topog- 
raphy of the region. The general circulation was studied using a two-lavered numerical 
model. Forcing of the model was done by varying: the inflow and outflow velocity; the 
port locations and boundary conditions; and by adjustment of topography. Represen- 
tative topography of the Fram Strait was modeled in the lower layer. Results indicate 
that the model is sensitive to EGC width and its proximity to the western boundary; and 
that in the absence of EGC flow, the Return Atlantic Current does not exist. Investi- 
gation of model sensitivity to WSC inflow vertical shear led to the conclusion that when 
the WSC inflow is 10 cm/sec, in the upper layer and 5 cm/sec in the lower layer, a 


boundary trapped component of the WSC 1s generated. 
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I. INTRODUCTION 


A. THE FRAM STRAIT REGION 

The area between East Greenland and Svalbard is considered to be of extreme 1m- 
portance since it provides the primary connection between the North Atlantic Ocean and 
the Arctic Ocean. The complex circulation and intense dynamical activity has for many 
years been of great interest, both from a military and scientific standpoint. Figure 1, 
from Paquette er al., (1985), shows the general circulation and topography of the region. 
The East Greenland Current (EGC) and the associated frontal jet (see figure 1), some- 
times referred to as the East Greenland Polar Front, moves cold low-salinity water from 
the Arctic basin to the south along the Greenland coast. The West Spitzbergen Current 
(WSC) (see figure 1) carries warm saline water to the north, west of Svalbard. According 
to Gascard e1 al., (1988), the West Spitzbergen Current generally flows toward the 
northwest along the Yermak Plateau (see figure 3) west of Svalbard and splits into se- 
veral branches in the vicinity of 78° - 79° N with a branch continuing to the northeast 
and the others turning west, one being associated with the Molloy Deep (see figure 1) 
area (surrounded bv the Spitzbergen Fracture Zone) on the south, and the other asso- 
ciated with the Hovgaard Fracture Zone (see figure 3). Much of the definition of the 
circulation in the region has come from observations using shipboard sensors, the Mar- 
ginal Ice Zone Experiment (MIZEX) data sets, drifting buoys, and satellite imagery. It 
is generally thought that the West Spitzbergen Current. the jet-lhke East Greenland 
Current, and the Return Atlantic Current (RAC) (see figure 1) provide the basis for the 
Fram Strait circulation. This circulation generally follows the topography, Foldvik ez 
al., (1988), Gascard er al.. (1988), Paquette ez al., (1985), Quadfasel ez al., (1987). 

Figure 2 is from Aagaard er al., (1988), and shows the actual bathymetry of the re- 
gion in greater detail. A current meter was positioned at FS-9B (see figure 2) and the 
velocities obtained from FS-9B were compared with model current velocities (in a later 


section). 


B. THE OBSERVED VELOCITY AND TRANSPORT WITHIN THE FRAM 
STRAIT REGION 

The East Greenland Current is considered to be a western boundary current flowing 
generally southward with its strongest flow over the continental slope. The East 


Greenland Current has been observed to be baroclinic in nature, with currents typically 





Figure 1. Bathymetry and Currents in the Fram Strait (from Paquette et al., 
(1985)): Depth labels are in hundreds of meters. The 5600 meter de- 


pression at 79° North is the Molloy Deep. 


computed at 20-30 cm/sec, 1n the upper layer (Foldvik ev al., (1988)). Velocities as high 
as | m/sec have been reported by Paquette ef a/., 1985. This current produces a char- 
acteristic baroclinic transport of approximately 3 Sverdrups (Sv) over the upper few 
hundred meters. 

There 1s a predominant cyclonic circulation within the Greenland Sea, this being 


caused at least in part by the wind stress curl and the guiding of flow along the topog- 


Figure 2. 





VET. 


10° W 3° 


Fram Strait Bathymetry (from Aagaard et al., (1988)): © Depth labels 


are in hundreds of meters. Note position of 7S-9B. The rectangle tndt- 


cates the numerical model domaim (322 km by 395 km) used tn this 
study. 


raphy of the region west of Svalbard (Quadfasel er al., (1987)). The northern portion 
of this cyclonic gvre forms the Return Atlantic Circulation (RAC). Current speeds of 
up to 40 cm'sec have been observed in the West Spitzbergen Current, with a mean 
transport value of about 3 Sv (varied from 0 to 9 Sv), Aagaard (1982), Hanzlick (1983) 
Gascard et al., (1988). Aagaard ez al., (1988) compiled current meter data from the Fram 
Strait, and data from their FS-9B current meter (see figure 2.) will be compared with 


model currents in a later section. 


C. PURPOSE OF THIS STUDY 

The objective of this study is to examine how barotropic and baroclinic current 
flows interact with bathymetry to produce the observed Fram Strait circulation. This is 
accomplished through the use of a nonlinear, two-layer, regional numerical model with 


currents over idealized Fram Strait topography. 


Il. NUMERICAL SIMULATION OF THE FRAM STRAIT REGION 


A. THE NUMERICAL MODEL 
1. Model Equations 

The circulation within the Fram Strait is modeled using a two layer, semi- 
implicit, primitive equation, numerical scheme. The scheme was initially used by 
Hurlburt (1974) in ocean circulation studies and has been emploved numerous times 
(Hurlburt and Thompson (1980, 1982); Smith and O’Brien (1983); and Smith and Davis 
(1989)). Linear test cases have been run for comparison with linear analytic solutions 
to show model validity (Smith and Reid, 1982). Motion within each layer is governed 


by a momentum equation (2.1) and by a continuity equation (2.2). 
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The EGC jet (when included as a initial condition) is defined: 
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Subscript i denotes upper (i= 1) or lower G= 2) layer. V, and 1, are depth integrated 
transports in each laver. Upper layer mean thickness (//,) is chosen to be 200 m. The 
lower laver mean thickness (/Z,) is 5400 m. The first internal Rossby radius of defor- 
mation (R,) sssociated with this laver thickness distribution is approximately 14 km. 
The Coriolis parameter / is taken to be a constant. A, and A, are interface and surface 
distortion amplitudes (= 75.0 meters and 0.2 meters respectivelv). JL, is the e-folding 
width scale for the jet, and y, is the eastern edge of the jet located 67 km from the west- 
ern boundary. The amplitudes 4,, Aswere chosen to give a maximum Jet Velocity (tae) 
of approximately 40 cm’s in the upper layer and 10 cms in the lower layer in all exper- 


uments. 


Variables and notation are defined in the appendix. The fluid in both lavers 1s 
assumed to be hydrostatic and Boussinesq, and the fluid density tn each laver is constant. 
The effects of winds, ice, tides, thermodynamics and thermohaline mixing are not in- 
cluded. 

2. Model Domain 

A rectangular region (521.7 km x 395.0 km) was divided into 4.7 km by 5.0 km 
rectangles to form the grid for the numerical model finite differencing. The grid was 
oriented to have the Molloy Deep and nearby Hovgaard fracture zone near the center, 
at 79.40° N, 1.53” E. 

Topography 1s included by applying a smoothed field of gridded bathymetry into 
the model for each corner of the 4.7 km by 5.0 km rectangles. Due to a model constraint 
that the laver interface cannot intersect the free surface or the topography, the 
Shallowest topography was 600 meters. Because of this constraint, shallow topography 
and nearshore processes could not be included within this model. Figure 3 shows the 
topography used within the model, while Figure 2 shows the actual Fram Strait topog- 
raphy. Note that the major features of the topography like the Molloy Deep, Yermak 
Plateau, East Greenland continental slope, and fracture zones are preserved. Small scale 
features, such as several seamounts in the Molloy Deep, which are thought to be im- 
portant for the circulation there (Bourke er al., 1987), are not resolved. 

3. Boundary Conditions 

Two separate and distinct boundary conditions were used in this model. For the 
east and west sides of the Fram Strait a no-slip boundary condition was used. For the 
open (north and south sides) portions of the Fram Strait a prescribed inflow and outflow 
boundary condition was set to represent the West Spitzbergen current and the East 
Greenland Current. The inflow velocity was ramped up to its maximum value in ap- 
proximately 5 to 6 days in cases which are boundary forced only. In cases in which an 
EGC jet extends across the domain in the initial condition, the EGC inflow « -ndition is 
maintained constant in time. The lateral distribution of the inflow is Gauss.an with an 
e-folding scale of 20 km. For the EGC, the Gaussian function 1s centered in the inflow 
port; for the WSC, the jet is centered on grid point 10. An open outflow radiation 
boundary condition (Camerlengo and O’Brien (1980)) was specified on the north and 
south boundaries anvwhere that did not have inflow. Inflow and outflow specifications 


are described in greater detail in the experiments section. 
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Figure 3. Modeled Fram Strait Topography: Depth contours are indicated in 
meters. Locations of inflow boundary forcing are specified in terms of 


grid points on the north and south boundaries. 


4. Preliminary Considerations 
For a linear, inviscid, two-layer fluid, fluid particles can be shown to conserve 


potential vorticity 











Wess: 
QO; = 
h; 
7 Gy, Ou, a 
Where ¢, 1s relative vorticity (= ea Fi ) in each layer and 4h, 1s instantaneous layer 


thickness. 

Figure 4 shows initial upper and lower layer potential vorticity for the EGC jet 
initialized case at dav zero. The jet extends along the whole western boundary. Con- 
sidering the lower layer potential vorticity first, it 1s obvrous that the major contributions 
to Q, are the depth variations associated with bathymetry. Contrast this figure with 
figure 3 (topography). Thus, if frictronal or nonlinear effects are small, fluid particles in 
the lower layer will follow contours of bathymetry. There may be different circulation 
patterns possible depending on where a fluid particle enters the domain. The Rossby 


number for the lower laver of the initial EGC: 


U> ] 


TL ~ (fx 20,000) ~ 
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indicates that nonlinear eflects are unlikely to overcome the strong topographic steering. 
Likewise. frictional effects are small in the momentum balance (equation number 2.1) 
by choice of B, Transient vorticitv-conserving wave motions, such as topographic 
Rossby waves, can however also exist in which cross isobath changes in vorticity asso- 
ciated with changes in A, are balanced by relative vorticity changes ¢,. 
In the upper laver, the jet’s relative Vorticity dominates the potential vorticity 
QO,. In contrast to the lower layer, the Rossby number for the upper layer 
uy 
TL = .]4 
indicates that nonlinear effects may not be negligible. 
Associated with lateral and vertical shears in the initial conditions, conditions 
for barotropic and baroclinic instability mav be met. It is a goal of this research to ex- 


amine topographically steered flows in the Fram Strait region for stable flows before the 





Ficure 4. Upper and Lower Layer Potential Vorticity: North is to the nght in 
both charts. Upper and lower layer contour intervals are 25 x 10° 


Seca 


more complex problem of unstable flows in the region can be examined. The time scale 


for these instabilities is likely longer than the duration of the experiments presented here. 
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I. EXPERIMENTS 
A. PURPOSE OF THE EXPERIMENTS 


The purpose of the following experiments 1s to understand the factors that contrib- 
ute to the general circulation of the Fram Strait. A Key objective of this study is to de- 
termine the effect topography has on the circulation, and this is reasonably well 
demonstrated with the model. Additionally, questions of baroclinic versus barotropic 
inflow forcing are explored as well as initialization of the East Greenland Current with 
a jet. 

Output from the model is in the form of charts of upper and lower layer velocities, 
and charts of upper and lower layer potential vorticity. All simulations were run for 10 
davs. This period was adequate for topographically steered flows to evolve, but did not 
allow the flow to evolve to the point of demonstrating mesoscale instabilities. To aid in 
comparing modeled output with observed data, the current was sampled at a specific 
location within the model domain and then compared with im situ current meter data 
from Aagaard er al., (1988). The model current 1s measured at grid position (36, 56). (see 
figure 3 for location) within the model domain so that model output velocity at this lo- 


cation can be compared with FS-9B velocity data. 


B. PRELIMINARY EXPERIMENTS 

Table 1 lists the preliminary experiments and parameters that were varied. Exper- 
iments No. 1 through 3 were EGC boundary forced (inflow allowed only at a specified 
location and outflow everywhere else on the boundary), and Experiments No. 4 through 
6 were EGC jet initialized (equations 2.3 and 2.4). All of the prelinunary experiments 
were run utilizing the modeled Fram Strait topography of figure 3. All velocity fields 
are displaved at 10 davs. Shears, such as 40 cmisec in the upper laver and 10 cm/sec in 
the lower layer will be indicated by a shorthand notation of 40'10 -m/sec, or as 40 10 
cm'sec in the upper and lower layers respectively. This form of notation will be used tn 
figure captions as well as text. 

Each numerical simulation was run with a variation of one parameter. Of interest 
was how the East Greenland Current (EGC) and the West Spitzbergen Current (WSC) 
interact with the region’s topography to produce the observed topographically steered 


ions. In this section the results of the preliminary experiments are described. 


i 


1. Experiment No. 1 (No Jet Initialization, WSC shear 10/10 cm/sec) 

In this simulation there was no jet imitialization of the EGC, and the WSC was 
initially barotropic with a velocitv of 10 cm/sec in each laver. While a portion of EGC 
evolves along the western boundary, a substantial portion is steered down the Lena 
trough. Return Atlantic Current flow exists via the forcing of the WSC and the joining 
of EGC water flowing down the Lena trough. Large velocities are seen in the Hovgaard 
Fracture Zone associated with convergent effects of topography. A weak along bound- 
ary current component of the EGC is observed in the upper layer. No northward along 
boundary flow of WSC is observed. An anticyclonic feature is seen over the Yermak 
plateau in both the upper and lower layers and will be discussed later in this section. 


The velocity fields are depicted in figure 5. 


Table 1. PRELIMINARY VARIATION OF MODEL PARAMETERS 
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2. Experiment No. 2 (No Jet Initialization, WSC shear 10/5 cm/sec) 

The WSC inflow vertical shear was changed to 10,5 cm/sec in the upper and 
lower lavers respectively for this simulation. RAC exists via a westward flowing branch 
of the WSC in both the upper and lower layers. The formation of along boundary flow 
in both the WSC. .d EGC regions 1s seen in the upper laver and is suggestive of an 
eastern branch of the WSC and an along isobath jet-like structure to the EGC. The 
Velocity fields are depicted in figure 6. 

3. Experiment No. 3 (No Jet Initialization, WSC shear 10/1 cm/sec) 

In this simulation, the WSC inflow shear was again increased. The along 

boundary component of the EGC 1s seen in the upper laver. There appears to be no 


coupling of a WSC component with the southward arctic flow from the Lena trough in 


Figure 
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Experiment No. 1 (No Jet Initialization, WSC shear 10/10 
cm/sec): The upper and lower figures represent upper and lower laver 
velocity fields respectively for Experiment No. 1. North ts to the right 
Mmeteln tere aNd edciinhe mepiesents a o22 chm OY 395 kin area. 
Contour intervals are 10 and 20 cm/sce for the upper and lower laver 


respectively. 
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Figure 6. 
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Experiment No. 2 (No Jet Initialization, WSC shear 10/5 cm/sec): The 
upper and lower figures represent the upper and lower layer velocity 
fields respectively for Experiment No. 2. North is to the right and cach 
figure represents a 522 km by 395 km area. Contour intervals are 5 and 


10 cni/sec for the upper and lower lavers respectively. 


the upper layer, but the WSC appears to connect with the RAC in the lower layer. The 
upper WSC has formed an along boundary northward flow pattern. The velocity fields 
are depicted in figure 7. 

4. Experiment No. 4 (Jet Initialization, WSC shear 10/10 cm/sec) 

This simulation was initialized with an EGC jet across the domain (along west- 
ern boundary). Its orientation and location relative to the East Greenland boundary is 
shown in figure 4 (upper). Initially, the jet’s upper layer velocity was 40 cm/s and lower 
layer velocity was 10 cm/s. The EGC inflow was boundary forced with 40/10 cm/s in the 
upper and lower layers respectively. The WSC was boundary forced with 10/10 cm/s in 
the upper and lower layers. On day 10, along boundary EGC exists in the upper laver. 
The RAC is not well organized in the upper layer in the region of the juncture with the 
EGC outflow, but does exist. A weak, narrow, along boundary component of WSC is 
seen in both the upper and lower layers. The velocity fields are depicted in figure 8. 

5. Experiment No. 5 (Jet Initialization, WSC shear 10/5 cm/sec) 

As in the previous experiment an EGC jet was initialized, with velocities of 40 
cm’s upper layer and 10 cms lower layer. Weak EGC along boundary flow exists in 
both the upper and lower layers, but appears more organized in the upper layer. WSC 
induced RAC exists. The westward flowing RAC 1s not well organized near its junction 
with the EGC in the upper laver, but the lower laver appears more organized. The along 
boundarv component of the WSC 1s well formed in the upper layer, but verv weak if at 
all in the lower laver. The velocity fields are depicted in figure 9. 

6. Experiment No. 6 (Jet Initialization, WSC shear 10/1 cm/sec) 

This simulation was initialized with an EGC jet across the domain in the upper 
and lower lavers. Little along boundary continuous deep EGC 1s seen with the exception 
of a small amount near the southern outflow region. Instead. cross isobath topographic 
Rossby wave variability extends along much of the slope. WSC induced RAC does not 
exist in the upper layer. The WSC component is totally boundary trapped i the upper 
layer. The lower layer shows RAC to the west. There is no along boundary component 
of WSC in the lower layer. The velocity fields are depicted in figure 10. 

Preliminary Experiments No. 1 through 6 all exhibit an anticyclonic vortex or 
vortices in the vicinity of the Yermak plateau. The Yermak plateau is the shallowest 
region within the model domain with a-mean depth of 600-1200 meters. The boundary 
forcing of the WSC inflow at the 2400 meters to 3200 meters isobath (southern bound- 
ary) causes along isobath flow toward the Yermak plateau. To exit the open boundary 


in shallow water on the Yermak Plateau, fluid columns are forced upslope. This flow 1s 


—y 
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Figure 7. Experiment No. 3 (No Jet Initialization, WSC shear 10/1 cm/sec): The 
upper and lower figures represent the upper and lower velocity fields re- 
spectively for Experiment No. 3. North is to the right in each figure and 
each figure represents a 522 km by 395 km area. Contour intervals are 


5 cmi/sec for both the upper and lower layers. 


Figure 8. 
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Experiment No. 4 (Jet Initialization, WSC shear 10/10 cm/sec): The 
upper and lower figures represent the upper and lower laver velocity 
ficlds respectively for Experiment No. 4. North is to the right in each 
feunesanidecach Maure represciits a S22 kineby 395 kin area. Contour 
intervals are 10 and 20 cm/sec for the upper and lower lavers respec- 


tively. 
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Figure 9, 
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Experiment No. 5 (Jet Initialization, WSC shear 10/3 cm/sec): The 
upper and lower figures represent the upper and lower laver velocity 
fields respectively for Experiment No. 5. North is to the right in each 
figure and each figure represents a 522 km by 395 km area. Contour 


intervals are 5 and 10 cm/sec for the upper and lower lavers respectively. 
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Figure 10. 
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Experiment No. 6 (Jet Initialization, WSC shear 10/1 cm/sec): The 
upper and lower figures represent the upper and lower Javer velocity 
fields respectively for Experiment No. 6. North is to the right im each 
figure and each figure represents a 522 km by 395 km area. Contour 


intervals are 5 cm/sec for both the upper and lower layers. 


prohibited from following deep isobaths around Yermak plateau by EGC deep flow 
down Lena trough in the opposite direction. This upslope flow causes vortex squashing 
over the Yermak plateau and the anticvclonic vortex or vortices is/are produced. 
Follow-on experiments will investigate methods of removing these vortices, since they 
are artificially induced by the model constraints. 

When Experiments No. 1 through 3 are compared with experiments 4 through 
6 it 1s seen that the model is not particularly sensitive to EGC boundary forcing versus 
EGC jet initialization with boundary forcing. The most noticeable difference is a broad 
jet (in the jet initialized experiments) versus a boundary trapped jet (in the boundary 
forced cases). In jet initialization cases, the lower layer is rapidly eroded by topographic 
waves. During the early stages of these simulations, these topographic waves naturally 
arise in the lower layer as a result of forcing with the EGC jet initialization and boundary 
forcing. These waves have a period of approximately 4 davs and wavelengths of ap- 
proximately 50 km. They propagate along the slope and out the open boundaries within 
10 davs. 

A comparison of experiments | through 3 and of experiments 4 through 6 leads 
to the notion that the existence of the WSC in two branches (one feeding RAC to the 
West, the other supplying along boundary flow) is dependent upon inflow vertical shear 
(specifically, 10°5 cm‘sec in the upper and lower layers respectively was used). No 
northward branch of WSC 1s seen in barotropic cases. 

During 1985-1986, Fram Strait current measurements were obtained. Current 
meter measurements from FS-9B were reported by Aagaard e? al., (1988) and are in- 
cluded here in table 2. FS-9B data is compared with model current meter data for a 
model domain location (36, 56) coincident with FS-9B, listed in Table 3. FS-9B is in the 
RAC region. Bv sampling the model current velocity in the same location, the ability 
of the model to correctly simulate RAC is tested. It should be noted that the data at 
FS-9B was obtained over a record length as indicated in table 2, and the miodel current 
velocities are the value on dav 10. This necessarily implies that direct comparison can- 
not be made. However, a qualitative comparison between the two data Sets suggests 
that model output velocities are about twice that of observed. For example, if Exper- 
iment No. 3’s upper laver current velocity of 22 cm/sec 1s compared with the value of 
FS-9B’s current meter at 107 meters, a value of 9.9 is seen. The direction of the current 
sensed at current meter FS-9B is westward. as is that of the model's current. Add- 
tionally, the model flow in this region is strongly barotropic despite baroclinic inflow 


conditions in some simulations. The convergent effect of topography creates enhanced 


lower layer flows in this region. Lower layer velocity at this location is not sensitive to 
the lower layer inflow of WSC indicating that much of the flow there (at least in these 
runs) originates from deep EGC flow down the Lena trough. A method to reduce the 
modeled current velocities may be to move the inflow region of the EGC to the west, 
more up on the slope and thereby reduce the amount of EGC flow down the Lena 
trough. 


Table 2. FRAM STRAIT CURRENT MEASUREMENTS (1985-1986) 
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Table 3. MODEL CURRENT DATA AT 10 DAYS 





C. FOLLOW-ON EXPERIMENTS 

Follow-on experiments were conducted to address the following issues: can more 
realistic current flow be obtained at the model position (36, 56)?; 1s the EGC responsible, 
at least in part, for the existence of RAC as observed?; and, is the miodel sensitive to 
inflow location of the EGC?. 

Table 4 lists the follow-on experiments and parameters that were varied. Exper- 
iments No. 7 through 9 address the Yermak vortex issue, while Experiment No. 10 is a 
flat bottom (no topography) case. Experiments No. 11 through 13 address EGC port 
variation and responsibilitv for RAC. Experiment No. 13 has modified topography and 
appears to have many of the features seen in observations. None of the follow-on ex- 
periments contained ject initialization. Inflow locations are given in grid point format 
and the specific locations may be scen in the representative Fram Strait modeled to- 


pography for each experiment. The flat bottom (no topography) experiment 1s included 
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Experiment No. 7: Fram Strait Modified Topography 


Figure 11. 
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as an ald in concluding that flows such as the RAC and Yermak anticyclone cannot exist 
in the absence of topography. 
J. Experiment No. 7 (Topography Modified) 

Experiment No. 7 1s boundary forced with a shear of 40/10 cm/sec in the EGC 
between 55-75, and in the WSC with a 10’5 cm/sec between 5-35. The bottom topog- 
raphy was modified near the eastern boundary such that there were no depths shallower 
than 1200 meters (see figure 11), thereby moving the Yermak plateau depth to a mini- 
mum of 1000 meters below the interface depth. The major difference here is that the 
large velocity Yermak anticyclone no longer occurs. Jn the upper laver there is an or- 
ganized boundary flow along the EGC, RAC exists and appears to be driven by a com- 
bination of EGC flow down the Lena trough and WSC flow that has split and has a 
Westward component and a boundary trapped component. In the lower laver there is 
no EGC boundary current. RAC exists, and WSC boundary current exists but to a 


lesser degree than in the upper laver. The velocity fields are depicted in figure 12. 


Table 4. FOLLOW-ON VARIATION OF MODEL PARAMETERS 
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2. Experiment No. 8 (Topography Modified) 

As an alternative way of changing the Yermak topography. the 1200 meter 
isobath was drawn to the northeast corner of the domain (see figure 13) allowing deep 
outflow all along the northern boundary. The EGC was boundary forced with 40, 10 
Cmpsec and the WSC was boundary forced with 10/5 cm/sec in the upper and lower 
layers respectively. In the upper layer, boundary trapped WSC exists. There is no 


Menmak anticyclone im the upper laver: RAC ‘exists. In the lower laver there is no 
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Figure 12. Experiment No. 7 (Topography Modi? ed): ‘The upper and lower fig- 
ures represent the upper and lower layer velocity fields respectively for 
Experiment No. 7. North is to the right in each figure and each figure 
represents a 522 km by 395 km area. Contour intervals are 5 cm/sec 


for both the upper and lower lavers. 
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boundary trapped EGC, there is RAC, and there is no boundary trapped WSC. AI- 
though there is no specific Yermak anticyclone, there is a more complex circulation in 
the east central boundary area of the domain. These motions have strong cross isobath 
components, indicative of topographic Rossby or continental shelf waves. The velocity 
fields are depicted in figure 14. 

3. Experiment No. 9 (Topography Modified) 

In this simulation, the topography was modified (see figure 15) in the area of the 
WSC inflow such that isobaths were brought parallel to the eastern boundary, thus al- 
lowing along isobath inflow, Foldvik (personal communication, 1989). The WSC was 
boundary forced with 10,5 cm/sec in the upper laver and lower layers respectively. The 
EGC was boundary forced with 40/10 cm/sec in the upper and lower lavers respectively. 
In the upper laver there ts a well organized along boundary EGC, there is RAC fed bv 
WSC and EGC flow down the Lena trough, and there is WSC boundarv flow. In the 
lower laver there is some EGC boundary flow near the exit region, and westward flow 
(RAC) in the south that appears to be driven by southward flow through the Lena 
trough. WSC boundary flow exists in both layers. Since the Yermak region was not 
modified. we see its anticyclone in both he upper and lower layers. The velocity fields 
are depicted in figure 16. 

4. Experiment No. 10 (Flat Bottom) 

In this experiment the topography was removed giving a uniform bottom laver 
thickness of 5400 meters with an upper layer thickness of 200 meters. The EGC was 
boundary forced with 40:10 cm'sec and the WSC was boundarv forced with 10.5 cm sec 
in the upper and lower layers respectively. As expected. along boundary flow (east and 
mest) aeveloped in the upper layer and lower layers. The RAC and the Yermak 
anticyclone are not formed in the absence of topography. The velocity fields are depicted 
imenecure 17. 

5. Experiment No. 11 (No EGC) 

In this simulation the EGC boundary forcing was eliminated. The WSC was 
boundary foreed with 10.5 cm/sec im the upper and los.er lavers respectively. The enure 
northern boundary is open to allow for outflow. In the upper laver there 1s northward 
EGC flow, reverse RAC flow and this couples with WSC to flow north through the Lena 
trough and out the northern boundary. In the lower laver there 1s sinular reverse flow. 
This demonstrates the need for the EGC to flow southward, thereby causing westward 
flow of RAC flow as seen in observations of Aagaard er a/., (1988). The velocitw fields 


are depicted in figure iS. 
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Figure 13. 9 Experiment No. 8: Fram Strait Modified Topography 
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Figure 14. Experiment No. 8 (Topography Modified): The upper cod lower fig- 
ures represent the upper and lower layer velocity ficlds respectively for 
Experiment No. 8. North is to the right mm each figure and each figure 
represents a 522 kin bee 395 kin area. Contour intervals are 5 cm, sec 


for both the upper and lower layers. 
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Figure 15. 9 Experiment No. 9: Fram Strait Modified Topography 
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Figure 16. Experiment No. 9 (Topography Modified): The upper and lower fig- 
ures represent the upper and lower laver velocity fields respectively for 
Experiment No. 9. North is to the right im each figure and each figure 
represents a 522 km by 395 km area. Contour intervals are 5 and 10 


cni’/sec for the upper and lower lavers respectively. 
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Fieure 17. 


Experiment No. 10 (Flat Bottom): The upper and lower figures rep- 
resent the upper and lower layer velocity fields respectively for Exper- 
iment No. 10. North is to the right in each figure and each figure 
represents a 522 km by 395 km area. Contour intervals are 5 and .5 


cm/sec for the upper and lower lavers respectively. 
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Figure 18. Experiment No. 11 (No EGC): The upper and lower figures represent 


the upper and lower layer velocity fields respectively for Experiment 
No. 11. North is to the right in cach figure and each figure represents 
a 522 km by 395 km area. Contour intervals are 5 and 10 cm/sec for 


the upper and lower layers respectively. 
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6. Experiment No. 12 (EGC Port Closer to Western Boundary) 

In this simulation the EGC was boundarv forced with 40/10 cm/sec at 65-75 in 
the upper and lower lavers respectively. The WSC was boundary forced with 10/5 
cm/sec at 5-35 in the upper and lower layers respectively. In the upper layer, EGC 
boundary flow exists, WSC induced RAC exists, and WSC boundarv flow exists. Little 
to no EGC flow is seen in the Lena trough. In the lower layer, WSC induced RAC is 
seen, and the beginning evolution of boundary trapped EGC is seen. The velocity fields 
are depicted in figure 19. 

7. Experiment No. 13 (Topography Modified, EGC port 65-75) 

This experiment was conducted with EGC boundary forcing of 40/10 cm/sec at 
65-75 in the upper and lower layers respectively. The WSC was boundary forced with 
10/5 cm'sec at 5-35 in the upper and lower layers respectively. The topography was 
modified (see figure 20) in two areas: first, the WSC inflow region in the south was 
modified to provide for along isobath flow of the WSC inflow (like that of Experiment 
No. 9); secondly, the Yermak plateau region was modified by drawing the 1200 meter 
isobath to the northeast corner as was done in Experiment No. 8. In the upper layer, 
EGC along boundary flow exists, RAC exists driven by a westward component of the 
WSC, and weak along boundary WSC is demonstrated. In the lower layer, along 
isobath flow is seen in the vicinitv of the EGC, RAC driven by WSC westward flow is 
seen, and a northward component of the WSC is seen just east of the Molloy Deep. 
Northeastward outflow is seen with slight WSC along boundary flow in the south. In 
the Lena trough, flow is southward on the western side and northward on the eastern 
side. The northward flow is joined from the south by a northward branch of the WSC 
that flows northward just east of the Molloy deep area. The velocity fields are depicted 
in figure 21. 

These flow fields are also seen in figure 22 which shows surface and interface 
height anomaly fields. A comparison of the surface slope field with the corresponding 
upper ocean velocity vectors indicates that the flow is geostrophically balanced, with 
flow to the right of the pressure gradient. Interface height anomalv ts indicative of the 
degree of baroclinicity. Largest interfacial slopes (and hence vertical shears) are seen in 
the EGC region where a downslope of 65 meters in the interface is seen toward the west. 
The RAC is relatively barotropic. 

Although Experiment No. 13 appears to best approximate observed flow di- 


rections, it is important to now consider magnitudes of velocities within its domain. For 
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Figure 19. Experiment No. 12 (EGC Port Closer to Western Boundary): The 
upper and lower figures represent the upper and lower laver velocity 
fields respectively for Experiment No. 12. North is to the right m each 
Hetinesana each (eunemepresents alone Kinkoy 395 Kil area, Contour 
intervals are S and JO cm/sec for the upper and lower lavers respec- 


tively: 
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experiment No. 13: Fram Strait Modified Topography 


Fieure 20. 
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Figure 21, 
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Experiment No. 13 (Topography Modified, EGC port 65-75): The up- 
per and lower figures represent the upper and lower Jayer velocity fields 
respectively for Expernnent No. 13. Northis to the nght in each figure 
and each hgure represents a 522 km by 395 km area. Contour intervals 


are 2.5 ci/sec for both the upper and lower layers. 





Pictine =. 2. 


Height Anomaly Fields: The upper field represents the surface height 
anomaly (Contour Interval = 2.5 cm). The lower ficld represents the 
interface height anomaly (Contour Interval = 10m). Positive (nega- 
tive) contours represent upward (downward) distortions in the inter- 


face. 
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comparison purposes, table 5 displays model current speeds for selected follow-on ex- 


periments versus current meter data from FS-9B presented in table 3. 
Table 5. MODEL CURRENT DATA AT 10 DAYS 
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Of particular note is that even for a qualitative comparison, Experiments No. 


12 and 13 have velocities close to those obtained at FS-9B, with Experiment No. 15 be- 
ing very close to the same. This coupled with the flow direction similarity of Experiment 
No. 13 is suggestive of a correctly mode!:d Return Atlantic Circulation. 

The follow-on experiments indicate a number of more realistic aspects of Fram 
Strait circulation than do the preliminary experiments. With respect to the RAC, there 
is reduced magnitude of the flow in Experiments No. 12 and 13 where EGC flow is 
forced only on the western side of the Lena trough. With respect to the Yermak 
anticyclone. its absence when shallow isobaths do not extend to the northern boundary 
(Experiments No. 7, 8 and 13) indicates an improper choice of model domain location 
relative to isobaths in early simulations. A more appropriate boundary location would 
follow isobaths. Irregular boundaries are not however allowed with the present model. 
A shift of the present domain 1° north would allow deep along isobath WSC inflow. and 
deep along isobath outflow (as seen in Experiment No. 13). In Experiment No. 13, the 
WSC can realistically split near the Molloy Deep into westward RAC with a small 
northward flow eas. of the Mollov Deep. A portion of this flow could then go around 
the Yermak plateau to the northeast, and a portion could exit to the north via the Lena 


trough. 


IV. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The Fram Strait with its associated East Greenland Current (EGC), West 
Spitzbergen Current (WSC), and varied topography ts a region of complex circulation 
and intense dynamical activity. To better understand the topographically steered flows 
in the Fram Strait will provide one with the ability to better and more effectively employ 
Naval forces throughout the region and to allow for continued scientific research. 

A two layer, semi-implicit, primitive equatron, numerical model is used to simulate 
the general circulation within the Fram Strait. The stmulation of the EGC and WSC 
using boundary forcing provides a more realistic circulation pattern. The model is not 
particularly sensitive to jet versus no jet initialization. Although the EGC ts observed 
to be jet-like (in the area of the East Greenland Polar Front) tn character, from a nu- 
merical modeling standpoint jet mitialization is not necessary for correct model simu- 
lation of the experiments of this study. 

The model ts sensitive to WSC tnflow vertical shear. In simulations with no shear, 
boundary trapped flows on the eastern and western boundaries could not be produced. 
When the WSC inflow vertical shear was increased to 10/5 cm/sec, the eastern and 
western boundary trapped flows were such to cause a reduction in model velocity, and 
thereby a tendency toward values more in line with those observed by Aagaard er ail, 
(1988) at FS-9B. Conversely, if shear 1s further increased such that upper layer flow 1s 
10 cm/sec and lower layer flow is 1 cm/sec, no appreciable reduction of model current 
Velocity is seen. Therefore, a shear of 10’5 cm/sec appears best to simulate RAC flows. 

Varving the port width of the EGC inflow allows for narrow along boundary flow 
of EGC or if widened a more broad flow extending south through the Lena trough re- 
gion. By moving the EGC inflow closer to the western boundary, there is a4 reduciem 
of flow down the Lena trough. This reduced flow down the Lena t ugh causes a re- 
duction of the model current velocities associated with the RAC region. EGC port 
forcing on the western side of Lena trough gives more realistic RAC flow values. Un- 
fortunately, no deep current meter measurements in Lena trough are available for com- 
parison with this model. 

The EGC 1s largely responsible for the westward RAC circulation in this model. 


Although some WSC ts required to substantially duplicate observed flow patterns within 


the general Fram Strait circulation, the major contributor is the EGC. This is demon- 
strated by the RAC reversal when EGC is absent. 


B. RECOMMENDATIONS 

A number of previous studies have indicated the importance of topography in 
Steering ocean flows in the Fram Strait region. This study illustrates the response of 
topographic steering to flow location and vertical shear of the currents. A number of 
other effects are however known to be important in the region. 

Further studies should incorporate the effects of winds and the marginal ice zone. 
The existence of an along ice edge wind driven jet could thus be included. Longer sim- 
ulations including this effect could address instability issues associated with the jet. The 
effects of mesoscale eddies on the general circulation could then be determined. 

In two-layer flows over topography, the upper layer may respond too strongly to 
topographic steering of the lower laver. This effect could be determined in future studies 
by the incorporation of a third model layer. 

In this study, inflow port locations were the same for upper and lower laver flows. 
The effect of forcing upper and lower flows at different locations should be addressed in 
future studies. Specifically, a broad upper laver East Greenland Current forced over a 


narrow deep inflow west of the Lena Trough may give more realistic results. 


V. APPENDIX- SYMBOLS AND NOTATION 


Biharmonic friction coefficient = 0X 10 es 

Initial surface and interfacial height anomaly of the jet 

Coriolis parameter for mean latitude 80°N 

Gravitational acceleration 

Reduced gravitational acceleration = g(p,— ,)/p, = .02 m/s? 
Instantaneous upper (1= 1) and lower (1= 2) layer thickness 


Upper (i= 1) and lower (1= 2) layer mean thickness 


Grid indices in x.¥ directions = 111,79 
e-folding scale for the jet = 20 km 
Pressure in the upper layer = g(, +h, +d) 
Pressure in the lower layer =p,-gh, 


Potential vorticity (f+ ¢,) / h, 

First internal Rossby radius of deformation 

= I/fe(e 7, AiG + 7) = 13.7 km 
Rossby number = 1 ax/ fL) 
Velocities in the x and v directions 

Transport in the x and y directions 

Maximum jet velocity 


Cartesian coordinates directed E and N respectively 


Veliax = 4.7 km 
Delta v = 3.0 km 
Model time step = 2400 sec 


Densitx ini 1? ler 
Gradient operator =6/ex+ e/éey 
Laplacian operator = 0/0 + 6 /ex 
Upper (1= 1), lower (1= 2) layer relative vorticity 
=¢y/¢x — éu/éy 


Eastern edge of the jet 
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